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Abstract 

We consider the structure of the leading ultra-violet (UV) renormalon 
singularity associated with the QCD vacuum polarization Adler D-function, 
in the approximation that only planar Feynman diagrams are retained. This 
“planar approximation” results in some simplihcations, in particular three 
of the four potential contributions from four-fermion operators are shown to 
be absent. We are able to obtain a fully normalized result for the leading 
n ^ oo behaviour of the portion of perturbative coefficients proportional to 
for SU(iV) QCD with Nf quark flavours. 
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The large-order behaviour of QCD perturbation series and the structure 
of the renormalon singularities which give rise to it have attracted a great 
deal of recent research interest [Q. For the QCD vacuum polarization Adler 
D-function the leading large-order behaviour is generated by an ultra-violet 
(UV) renormalon. The structure of this Borel plane singularity can be identi¬ 
fied following the work of Parisi p|, by considering the insertion of dimension- 
six operators into the Green function which generates D [^. In this paper 
we wish to make use of the results of ref. [3] to study the asymptotics of the 
perturbative coefficients of D in the “planar approximation” where at each 
order in the 1/A^/ expansion in SU(A^) QCD with Nf quark flavours, only 
the leading contribution in the large-A^ limit is retained. This corresponds 
to the evaluation of planar Feynman diagrams |^, and in this approximation 
some simplifications occur. 

The Adler D-function is defined in terms of the correlator 11 (s) of two 
vector currents in the Euclidean region, 

^{s){qt,qu - giMuq"^) = J d^x e*‘^'^(0|T{ Ji.(0)}|0) , (1) 

where s = —q^ > 0. By taking a derivative in s one avoids an unspecified 
constant and defines 

D(s) = -s^n(s). (2) 

In SU(A^) QCD perturbation theory, 

= + + b, (3) 

where Qf denotes the electric charge of the quarks and the summation is 
over the flavours accessible at a given energy. The SU(A^) Casimirs are 

Ca = N, Cf = {N'^ — l)/2N. D denotes the contribution of “light-by-light” 
diagrams which first enter at 0(af) in perturbation theory. We shall be 
interested in the corrections to the parton model result, D, which in “planar 
approximation” have the form, 

D = y; + O , (4) 
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where do = 1 and a = ag/n is the RG-improved coupling. We stress that 
the dn correspond to the perturbative coefficients in “planar approximation”. 
Crucially each dn is then a sum of multinomials in iVj and N of degree n, 


dn = ^ 

n 

= Y. , (5) 

r=0 

where the d]^~'^\ d^n~^^ are pure numbers, b = (lliV — 2iV/)/6 is the hrst 
QCD beta-function coefficient. The second expression is the “^-expansion” 
[5-7] obtained by rewriting Nf = (lliV — 6&)/2. 

Our aim is to obtain asymptotic results for and for hxed r and 
large n. Representing Z) as a Borel integral, 

D = r dze-^/'^Bimz) . ( 6 ) 

Jo 

The Borel transform B[D]{z) is expected to have a leading UV renormalon 
singularity aX z = —2/6. As noted earlier this should receive contributions 
from the insertion of dimension-six operators, so that 


, Ki -|- 0(1 -|- 6 ^/ 2 ) 

^ ^ i (l + 6^/2)'* ’ 

where the sum is over the contribution of operators Oi 
have the form 


B\D]( 


( 7 ) 


The exponents 6i 


Si=pi- (2c/6) -h Aj , 


( 8 ) 


where pi is an integer, c denotes the second universal beta-function coefficient. 


c = 


7^ 
8 b 


11 CaCf 
8 b 


-Ca + -Cf 


(9) 


and Xi denotes the one-loop anomalous dimension of the operator. 

In the “planar approximation” the residues Ki and exponents 5i can only 
involve the ratio N/b in order to reproduce the multinomials b'^~'^N'' in eq.(5). 
That is we must have. 
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where the and are pure numbers. The combination A ^/6 is the 

organizing parameter of an expansion in the number of renormalon chains 
exchanged in QCD planar diagrams, and replaces the 1/iV/ appropriate to 
the chain expansion in QED [^. 

To obtain the large-iVj asymptotics we need to consider the one-chain 
contribution generated by the operator Oi = {1/g^)dijF'^^^dPFp^^ with g the 
QCD coupling and Fp_^ the held strength tensor associated with the external 
vector source. The large-A^/ behaviour is of course known exactly i 0 . 
The one-loop anomalous dimension of Oi vanishes and 


5i = 2- ( 2 c/ 6 ) . 


( 11 ) 


If we work for convenience in the “V-scheme” (MS with scale fF = e 
then one hnds /T® = 4/9 Si • Exchanges of two and more renormalon 
chains will involve the four-fermion operators Oy and Oa P, , 


Oa = (V’7m75^^^)(V’7^752"^^/’) , (12) 


where denotes colour matrices (these correspond to 03,04 in ref.[3]). 
The remaining four-fermion operators which do not contain ( 0 i ,02 in 
ref. [3]) will not contribute in the “planar approximation” since they will be 
suppressed by powers of N. 

The corresponding anomalous dimensions of the diagonal operator com¬ 
binations of Oy and Oa can be obtained from the anomalous dimension 
matrix computed in ref. [3]. On discarding contributions which vanish in 
“planar approximation” one hnds. 


Nj ^16iV2/9 + 9Ar2 _ 3JY 
“ ^ 46 


(13) 
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In the large-iV limit the diagonal operator combinations become 0± = 
Ov^Oa ■ After Fierz-ing O- is purely scalar and will not contribute to the 
vector correlator. This can be seen by showing that the large-behaviour 
involves traces containing an odd number of gamma matrices , which there¬ 
fore vanish. Thus = 0, m = 1, 2, 3,..., and there is no branch point 
with exponent in “planar approximation”. 

Thus in “planar approximation” the relevant Borel plane singularities are 
the one-chain contribution with branch point exponent and the 0+ con¬ 
tribution with exponent 6+. From the explicit large-two-chain calculation 
of [1^ one can deduce = —1/18 (in the V-scheme). The computation of 
the full residue iF+ would entail a summation of all orders in N/h , requiring 
exchange of any number of chains in planar diagrams to be calculated. 

The large-order behaviour resulting from the contribution in eq.(7) is 


r(5.) 



^? 7,!(1 -|- 0{l/n)) . 


(14) 


The asymptotics are then determined by the operator Oi with largest <5*. 
In both the full theory and in “planar approximation” Oj^ is the dominant 
operator. One has 0 


= 3 - (2c/6) + A± , (15) 

where A+ > 0 and A_ < 0 for A^j > 0, 6 > 0, (A+,A_ correspond to Ai, A 4 
tabulated in ref. [3] for various values of Nf and N = 3). decouples in 
“planar approximation” as we have seen, and gives only a sub-leading contri¬ 
bution in the full theory. Evidently from eq.(ll) the one-chain contribution 
is subleading by at least 0{l/n). Other held operators are also sub-leading 
in the full theory 0 . 

Whilst the asymptotics of the full coefficient dn has an undetermined 
normalization we shall now show that one can obtain a fully normalized result 
for the asymptotics of the coefficients in “planar approximation”. To 

derive this we begin by considering the large-A^/ limit . In this limit A+ = —2 
from eq.(13) , and Ai = 0, so that the one-chain contribution dominates by 
0 (n). 

Expressing 5i in eq.(ll) as an expansion in N/b using eq.(9) for c one 
hnds 




13 N 25 


(16) 
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Setting = 4/9 and systematically expanding ^/r(5i) in eq.(14) in 
powers of N/h one then obtains 


dri 


nn\ 


^ ^ “ T*'“" “ I + 


(17) 


where 7 = 0.57722... is the Euler-Mascheroni constant which arises in the 
expansion of r((5i). The asymptotics of the sub-leading “ 6 -expansion” coef¬ 
ficients can be obtained from successive terms in the square bracket 

in eq.(17). Their evaluation would require the calculation of the for 

m < r, that is the insertion of the operator Oi into planar diagrams with at 
least r -|- 1 renormalon chains. 

The large-order behaviour of the coefficients on the other hand is 

determined by the “leading- 6 ” 6 ” term in eq.(17). Expanding this using the 
binomial theorem gives a contribution n'^N^NJ~^/r\, which will dominate by 
at least 0(n) over the contributions arising from the (iV/ 6 )”^ terms in the 
square bracket. Thus one obtains the fully normalized asymptotic behaviour 
(in the V-scheme) 




-n 


(1 + 0 


n 


(18) 


,+ird 

9' r! V 2 

Notice that for eq.(18) to hold it is crucial that, as argued above, there 
is no 0_ contribution in “planar approximation”. In the large-iVj limit 
A_ = 0 and since 0_ hrst potentially contributes at the two-chain level , 
and is thus 0{1/Nf) suppressed relative to the one-chain contribution, both 
should contribute equally to the asymptotics of the sub-leading , r > 0 . 

Thus if were not to vanish, the 4/9 in eq.(18) would be replaced by 

I H—^ j. The result in eq.(18) was derived in ref.[5] assuming only the 

one-chain contribution to the leading UV renormalon singularity. We have 
shown here that it remains valid on including any number of chains. 

We now consider the large-iV limit and the asymptotics of d|[l In this 
limit (9+ is the dominant operator and A+ = 0 from eq.(13). Thus the one- 
chain contribution is sub-leading by 0(l/n). Returning to eq.(5) we note 
that it can be rewritten in yet a third form. 


djj 


00 

6”-”iV;:d<"-”> 

r=0 



(19) 
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so we can obtain an alternative “6-expansion” by writing N = {6b + 2Nf)/ll. 
Correspondingly we can write the residnes Ki and exponents Si in eq.(7) as 
expansions in Nf/b. 



where the hat distinguishes these coefficients from the N/b expansion in 
eq.(lO). Whilst the N/b expansion corresponds directly to evaluating planar 
diagrams with an increasing number of renormalon chains, the Nf/b expan¬ 
sion has no such direct diagrammatic interpretation. Evaluation of 
requires an all-orders resummation of the N/b expansion. For instance 

= (21) 
m=l VJ-J./ 

As noted earlier from the exact large-iV two-chain calculation of ref. [12] one 
can infer = —1/18, so with our present knowledge iC® = (—1/33) -|-... 
, it is of course not even guaranteed that the sum in eq.(21) converges. 

On expanding eq.(13),and eq.(9) for c, in Nf/b one finds 

, 261 263 Nf 3680 TV? , , 

“ m + tST + + ■ ■ ■ ■ 

Inserting these Nf/b expansions into eq.(14) gives, analogous to eq.(17). 


dri 


K 


( 0 ), 


-b/2fn 


n^l40/121^j 


r( 261 / 121 ) 


1 + 


Nf 263 ,, ,. 261 ,,, 

H-(Inn - M -)) -h ... 

6 ViT® 726^ 


(23) 


where i/{x) is the logarithmic derivative of r(a;) which arises in the expansion 
of r((5_|_) . In principle the asymptotic behaviour of the coefficients in 

the alternative “6-expansion” in eq.(19) can be derived from the successive 
terms in the square bracket in eq.(23). However, as pointed out above, the 
evaluation of the requires a resummation of any number of renormalon 
chains. 
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The large-order behaviour of the coefficients will again be determined 
by the “leading-6” 6” term in eq.(23). On expansion using the binomial 
theorem one obtains a contribution , which will dominate by 

at least 0(n) over the {Nf/h)^ terms in the square bracket. So that the 
leading large-order behaviour is 



K. 


( 0 ) 


r (261/121) 


n 


r+140/121 


n\ 


t\ 


llx*^ 

12j 



(l + O 



( 24 ) 


Whilst the evaluation of would require the all-chains resummation 
of eq.(21) by taking ratios we can arrive at the fully normalized result 


4' 




(l + O 


Similarly from eq.(18) for one can obtain 



(25) 


jk-s] 

tin 




(i + o 



(26) 


The results in eqs.(25),(26) hold universally for any QCD observable 
whose large-iV and large-iVj asymptotics are each dominated by a single 
branch-point (operator) in the Borel plane. For instance in the present 
case of the vector correlator we have seen that the large-iV asymptotics are 
dominated by 0+ ,and the large-iVj asymptotics by the one-chain contribu¬ 
tion. This ensures the weak dominance of the large-iV and large-iVj- derived 
“leading-6” terms evident in eqs.(18) and (24), that is on binomial expansion 
they generate the correct leading asymptotics for and d|/“”l 

An amusing speculation is that if one were able to perform the all-chains 
summation in eq.(21) one would hnd = 0, then assuming that the 

0(1 + hz/2) term in the numerator of the 6+ contribution also vanished in 
the large-iV limit, the leading large-iV asymptotics would in fact come from 
the one-chain contribution. If the two and higher chain contributions to the 
residue Ki also summed to zero in the large-iV limit. 




(m) 


m=l 


6 

n 


m 


0 , 


( 27 ) 



then the standard large-iVj derived “leading-6” term ,which is clearly exact in 
the large-iV/ limit , would also give the correct leading asymptotic behaviour 
in the large-iV limit , up to effects dependent on the second beta- 

function coefficient which vanish in large-iVj but contribute in large-iV. This 
would give the “leading-6” approximation [5-7], which has found extensive 
phenomenological application [5,7,13-15], a much stronger motivation. 

We stress that this is pure speculation. The only fact at our disposal is 
the two-chain approximation to derived from the calculation of ref. [12], 
^ —1/33. This small and negative value makes it not inconceivable that 
77_|_ has a zero at iVj = 0 (equivalent to the large-iV limit iV/6 = 6/llin “pla¬ 
nar approximation”) and is thereafter positive for Nf > 0 in asymptotically- 
free QCD. 
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